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Abstract 

We use a recently developed lattice model to study the percolation of particles of different sizes 
and shapes in the presence of a polymer matrix. The polymer is modeled as an infinitely long chain 
to simplify the calculation but we make it more realistic by considering it semi- flexible. We study 
the effects of the stiffness of the polymer, the size disparity of the filler particles, their aspect ratio 
and the interactions between fillers and polymer on the percolation properties of the system. The 
lattice model is solved exactly on a recursive square Husimi lattice, which is an approximation for 
a square lattice. The solution represent an approximate solution for a square lattice. Our results 
are able to reproduce most of the experimental findings that have been observed in the literature. 
In particular, we observe how an increase in the size disparity of the filler particles dispersed in 
the matrix as well as an increase of their aspect ratio decreases the percolation threshold. 
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I. INTRODUCTION 



In the past several years, many studies have been devoted to the analysis of the physical 
properties of filled polymers. As explained in the previous papers , polymers are filled 
because the presence of fillers can improve various physical properties: mechanical and 
rheological properties can be modified, as well as the electrical resistivity, the absorption 
of ultraviolet and infrared radiation, the optical clarity, etc. Therefore, the fillers that are 
introduced in the matrix of a polymer composite serve a useful purpose. 

Mechanical properties have been studied more extensively. For example, the addition of 
different volume fractions of carbon black to a polypropylene (PP) matrix can enhance the 
Young's modulus of the matrix. The composite that is produced is consequently stiffer than 
the matrix j^|. The processing conditions can also affect the properties of a composite. For a 
fixed fraction of black filler, samples prepared by injection molding have a Young's modulus 
that is consistently higher than the modulus of samples prepared with compression molding. 
This effect has been attributed ^ to a better distribution of the filler in the PP matrix due 
to additional compounding under shear in the injection-molding equipment. Orientation 
due to elongational flow during injection molding might also play a role while compression 
molded samples are supposedly isotropic 

The addition of a filler never improves all the properties of a matrix. Usually a decrease in 
the elongation at break of polymer samples filled with carbon black accompanies an increase 
in their Young's modulus. The composites that have a high filler content have a very low 
deformation at break. They break before they are able to achieve a significant degree of 
orientational reinforcement due to drawing during the tensile tests. The elongation at break 
of the compression molded samples, as well as their tensile strength jjj, is lower than that 
of the injection molded samples indicating the presence of larger or more frequent defects, 
like voids or cracks, as a consequence of a less perfect filler dispersion. The addition of a 
filler to a matrix can improve other properties. If the electrical properties of the same PP 
composites filled with carbon black are studied one observes a sharp increase in the electrical 
conductivity and a sharp decrease in the elongation at break, at about the same amount of 
filler, around 7% in weight P, 0] . 

The increase of the electrical conductivity and the decrease of the elongation at break 
are related, and may be associated with the formation of a more or less percolating filler 



2 



phase in the composite; this filler phase allows for a much higher electrical conductivity on 
the one hand, but also leads to the formation of more frequent failure sites on the other 
hand. The value of filler content at which the physical properties of a filled composite 
change dramatically is considered to be the percolation threshold. In this case, the value 
of the percolation threshold for the two processes is the same. However, if the changes in 
two different physical properties are associated with different microscopic mechanisms, the 
threshold associated with different physical properties should not coincide, as explained in 
the previous paper 3|. 

The results for the conductivity and the elongation at break of the same composites 
that have been prepared by injection molding show a much higher percolation threshold, 
about 16% in weight jjj. This observation is consistent with the idea that because of the 
different preparation conditions, the carbon black is less perfectly distributed in the PP 
matrix in the compression molded samples than in the injection molded samples. This 
allows for the formation of conducting channels formed by carbon black particles at much 
lower concentrations. As percolation is expected to be affected by size variations of filler 
particle, it is important to quantitatively assess their effects on percolation. The particles' 
size, shape and polydispersity after the composite has been processed matters most because 
these properties can be greatly affected by the processing itself jsj. Our results here will show 
that an increased polydispersity of the filler particles may be responsible for the reduction 
of the percolation threshold, and has been observed also in other works 

Although polymers are good electric insulators, they are usually poor thermal conductors 
and their dielectric properties are strongly affected by the presence of moisture. One solution 
to this problem is to make polymer composites with increased thermal conductivity. This 
has been achieved by adding proper fillers, characterized by a high thermal conductivity, to 
the polymer matrix. The thermal conductivity increases as the volume fraction of the filler 
increases in all predicted by many theories and observed in experimental studies 

. Carbon fibers increase the thermal conductivity of the polyurethane matrix 
more than the alumina particles. This is certainly related to the intrinsic properties of the 
reinforcing particles, the thermal conductivity of carbon black being higher than that of 
alumina. Also, both an increase in the size of the alumina particles and a decrease in the 
aspect ratio of the carbon fibers decreases the thermal conductivity. From the point of view 
of percolation, the decreased thermal conductivity should correspond to an increase in the 
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value of the volume fraction occupied by the filler at the percolation threshold. We have 
here another experimental evidence of how the size of the particles can affect the percolation 
threshold of a composite system. In general, an increase of the size of the particles increases 
the volume fraction of the filler at the percolation threshold. A decrease of the aspect ratio 
of the filler particles has a similar effect. 

The shape of the particles has also a critical effect on the electrical properties of the 
composite. Higher aspect ratios facilitate the formation of clusters and therefore lower the 
threshold concentration . The dielectric constant of composites samples containing asym- 
metric fillers usually exhibit a strong dependence on the frequency, increasing significantly 
as the frequency decreases. This significant enhancement of the dielectric constant values at 
low frequency can also be explained by the formation of filler clusters in the composites P| . 
At low frequencies, the surface charges that contribute to the electrical polarization are given 
more time to rearrange themselves on the clusters that have larger dimensions, generating 
greater polarizations. It is important to observe that in many experimental works, a change 
in particle shape can be associated with an increase in the polydispersity of the particles 
so that the observed increase in the dielectric constant, and consequently the decrease in 
the percolation threshold, could be ascribed to an increase in the polydispersity of the filler 
particles, as well as to a combination of the two effects. The addition of 20% volume of 
nickel to a polyethylene (PE) matrix affects the dielectric properties of the polymer matrix 
0. As expected, the dielectric constant values increase as the filler is added to the poly- 
mer matrix. The increase in the dielectric constant in the presence of metallic fillers has 
been explained as due to electrical polarization of the charges induced on the surfaces of the 
metallic filler particles ^]. The dielectric constant of composites with filaments and flake-like 
fillers is higher than the one of the composites filled with nickel powder and increases with 
increasing aspect ratio. The enhancement of the dielectric constant at higher aspect ratios 
las been explained as due to an increased easiness of formation of clusters of filler particles 
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The presence of filler particles, and their nature, can also affect the rheological properties 
of polymers. In general, the addition of a filler increases the viscosity of the polymer matrix 
[12I [13L Q]. Such increase is usually proportional to the fraction of filler present in the 
system until a critical value is reached where a sharp jump corresponding to the percolation 
threshold occurs. The limiting shear viscosity, as well as the yield stress, of carbon black filled 
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ethylene propylene diene monomer (EPDM) system composites depends on the fraction of 
black. The sudden increase in the shear viscosity around a value of 15phr (parts per hundred 
rubber) of carbon black, is explained in terms of the development of a trend in filler particle 
networking that becomes very prominent at what is thought to be the percolation threshold 



15 



16 



171 ] . The incorporation of fillers also enhances the low shear rate melt viscosity of 
EPDM. The narrowing of the difference in the shear viscosity values between compounds 
of different filler loadings at high shear rates occurs as a result of a progressive breakdown 
of the filler network structure 18[. As mentioned above, it is very complicated to measure 
and model the viscoelastic properties of composite systems since these properties depend 
on the dynamics of the system. The shearing rate at which the viscoelastic properties are 
measured, for example, is extremely important since high shear rate tend to largely nullify 
the effects of filler networking and the effect of filler-polymer physical interactions jl^j ]. In 
general, decreasing the size and the aspect ratio of fillers increases the viscosity jl3j ]. 

There is enough experimental evidence to conclude that the properties of the composite 
are affected by a large number of parameters: the size, shape and distribution of the rein- 
forcing particles as well as the interactions between the particles and the polymeric matrix. 
It is often difficult, if not impossible, to separate the effects due to the different parameters. 
To summarize the experimental results, we can say that the following general trends have 
been observed: 

• the percolation threshold appears to decrease as the size disparity between different 
particles present in the system increases, 

• the percolation threshold appears to decrease with the aspect ratio of the particles 
present in the system, 

• the nature of the interactions between fillers and matrix affects the percolation prop- 
erties of the system. 
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II. PERCOLATION OF SINGLE-SITE FILLER PARTICLES IN THE PRESENCE 
OF A POLYMER MATRIX 



A. Model 



In order to describe the percolation of particles of different sizes in the presence of a 



polymer matrix, we have used a model for a semiflexible polymer due to F 



has been extended and used previously in our research group 
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ory Il9|,|2fj, which 



2J]. This polymer 



model, originally used to describe an infinite molecular weight incompressible polymer, has 



been recently implemented to take into account the finite 
presence of free volume j^, l^j]. The model developed in 



eng th of many polymers and the 
22j is defined on a square lattice 



and solved exactly on a square Husimi lattice, part of which is shown in Figure d The 
Husimi lattice is built in a recursive fashion starting from a central plaquette. The different 
size of the plaquettes that belong to different generations is just a matter of convention: 
one might consider all plaquettes to have the same size but then a lattice would become 
impossible to draw because of the crowding that one has as one moves away from the origin 
of the lattice. 

In the present research, we extend this model. We consider a system made of a semiflexible 
polymer matrix, in the limit of infinitely long polymer chains, and filler particles of different 
sizes embedded in such a polymer matrix. Each monomer belonging to the polymers occupies 
a lattice site, and the lattice bond connecting two occupied nearest-neighbor sites represents 
a polymer bond in the polymer. The number density <pp of infinitely long polymer chains is 
kept equal to zero. In the following, we will often loosely speak of a single polymer, which 
also corresponds to <pp = 0, to describe this limit, even though the limit only ensures that 
the number of polymers does not grow as rapidly as the system size N, where N denotes the 
number of lattice sites. In this limit, either there exist no polymers if their monomer density 
is zero or they exist if their monomer density is non-zero. It is possible in the latter case to 
have polymers cover a finite non-zero fraction of the entire lattice. This is the reason why 
we can have configurations in which two polymer bonds are parallel to each other within a 
plaquette without being connected to each other, see below. As explained elsewhere j^j , the 
limit of infinitely long polymers (<f)p = 0) corresponds to a situation in which the end-points, 
whose density (= 2<pp) is also zero, do not have any effect on the thermodynamics of the 
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system. As a consequence one might think of this limit as useful for chains that are long 
enough so that the end points represent a negligible fraction of the polymeric monomers 
present in the system. 

Keeping (ftp = fixed throughout the study (see below to see how we ensure (ftp = 0), 
we have studied the percolation properties of systems containing up to two different kinds 
of filler particles. First, we analyzed the percolation of single-site particles which occupy 
one single site each on the lattice, just as a polymeric monomer. In the following, the 
single-site filler particles are going to be called solvent. To simplify the lattice model, we 
take the solvent and the polymeric monomers to have the same size. This restriction can 
be easily removed at the expense of making the recursive calculation more cumbersome. 
This will be done elsewhere. Then, in the following sections, we study the percolation of 
polydisperse fillers made of single-site solvent particles and square-shaped particles in one 
case and single-site solvent particles and star-shaped particles in the other. 

The original model used in the literature to describe the thermodynamics of semiflexib 
polymer chains was developed originally independently by Flory and Huggins 

This model was the basis for the Gibbs-DiMarzio theory of the glass transition that has 
represented for many years one of the most successful theories used to describe the glass 



transition in polymers 



271 . |28| . Unfortunately the approximations made in this theory 
have been proved to be qualitatively wrong |29|, |3C| |3l| , even though the resulting physics 
behind the entropy crisis seems to be valid not only for longpolymers but also for smaller 



molecules, as demonstrated recently by our group 
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22, l£2|. This has been one of the 



driving forces towards developing a new theory. Our current model is an extension of Flory's 
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model and his model is recovered in a particular limit 

In our model, a polymer chain is assumed to consist of n equal monomers, each of the 
same size as the solvent molecule, which we always designate as the smallest among all the 
filler particles which are present in the system. For simplicity (see above), we consider the 
solvent to occupy a single lattice site. This condition can be easily relaxed, but we will not 
do that here. Any larger filler particle is assumed to contain n' > 1 contiguous monomers so 
that it covers n! contiguous sites on the lattice. Each site of the lattice is occupied by either 
a polymer chain or a filler particle, and the excluded volume effects are taken into account 
by requiring a site to be occupied only once. The polymer chain occupies a contiguous 
sequence of all the lattice sites connected by polymer bonds, just as is the case of a filler 
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particle with n' > 1, without forming any branches or loops. 

For concreteness and ease of discussion, one can think of the system as if it were defined 
on a square lattice. At every site which is occupied by a monomer, the polymer chain 
can assume either a trans conformation (the conformation is related to the state of two 
consecutive bonds, and, consequently, of three consecutive monomers, when the consecutive 
bonds are collinear, or one of the two possible gauche conformations, when the polymer chain 
bends. For a semiflexible polymer chain, every gauche conformation has an energy penalty 
e compared to a trans conformation. This is the only energy present in the original model 
developed by Flory. We set the energy for a trans conformation to be zero. The inclusion of 
this single energy for a gauche bond is not satisfactory since it leads to a continuous melting 
transition that does not occur in nature 2^. In order to obtain a meaningful description of 
the physics of semiflexible polymer chains, it is necessary to include other conformational 
energies that were not included in Flory's model. We introduce two such energies: an 
interaction energy e p associated with each parallel pair of neighboring bonds, and an energy 
£h for each hairpin turn within each square of the Husimi lattice |22j. Corresponding to 
each quantity N g , N p and iVh, the number of gauche conformations, parallel bond pairs and 
hairpin turns, respectively, there is an independent activity w, w p and w^, respectively. The 
activities w, w p and are defined as follows: 

w = exp (—(3e) , w p = exp (— (3e p ) = w a , = exp(— fish) = w . (1) 

Here, j3 is the inverse temperature T in the units of the Boltzmann constant, and a = e p /e, 
b = Eh/e. The original model developed by Flory is obtained when the last two interactions 
are absent, so that w p = u\ = 1. It should be stressed that is the additional energy 
associated with the conformation, once the energy of the two bends and a pair of parallel 
bonds have been subtracted out. Both e p and are associated with four-site interactions, 
since it is necessary to determine the state of four adjacent sites to determine if a pair of 
parallel bonds or a hairpin turn is present. It has been seen that a and b represent possible 
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measures of the stiffness of the polymer chain 

In the presence of single-site solvent particles, the total partition function of the system 
can be written as 

Z N = ^^^ c V^«f h > (2) 
where the number of solvent molecules N$ is controlled by the solvent activity r] = w~^, 



where /x is the chemical potential of the filler particle. The number P of polymers is con- 
trolled by another activity given by H 2 . The interaction between nearest neighbor pairs 
N c of filler particles and monomers of the polymer is given by the Boltzmann weight 
w c = exp (—/3e c ) = w c . In all the cases studied during this research, the polymers have 
been considered of infinite length, so that H = 0. 

B. States of a site 

As the model cannot be solved exactly on any regular lattice, we replace the original 
square lattice with a Husimi lattice [lj]. We consider a linear polymer that covers all the 
sites of the Husimi lattice along with solvent molecules which occupy one site each. In 
order to be able to describe all the phases present in the phase diagram of the system, it is 
necessary to use a labeling scheme in which, if the base site of the plaquette has the label 
m, the two middle sites have the label (m + 1) and the peak site has the label (m + 2). 
This is a slightly different labeling than the original one used in the previous paper, where 
all the sites, other than the base one, were labeled as (m + 1) sites. As explained in [jj, the 
difference between the two labeling schemes is related to the different scheme adopted for 
the solution of the equations that we need to solve in order to obtain the thermodynamics 
of the system. In this case we are going to use a two-cycle recursion scheme so that the 
current labeling becomes fundamental to properly describe the solutions. 

Consider an mth level square. The base site can assume five different possible states 
depending on the state of its lower bonds connected to the base site, as shown in Figure |2J 

1. In the I state, both lower bonds are occupied by the polymer chain. Since the polymer 
is linear, the two upper bonds in the lower (m — l)th level square must be unoccupied 
by the polymer. 

2. In the O state, both lower bonds are unoccupied by the polymer chain but the upper 
bonds in the (m — l)th level square are occupied so that a monomer belonging to a 
polymer chain occupies the base site of the plaquette under investigation. 

3. In the L state, only one of the lower bonds is occupied and the polymer occupies the 
left upper bond of the lower (m — l)th level square (we always think about left and 
right moving down the lattice towards its center). 
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4. In the R state, only one of the lower bonds is occupied and the polymer occupies the 
right upper bond of the lower (m — l)th level square. 

5. In the S state, the base site of the plaquette under investigation is occupied by a 
solvent particle. 

For m — 0, the lower square in the above classification is the square on the other side of 
the origin. 

We are interested jl| in the contribution of the portion of the mth branch C m of the lattice 
to the total partition function of the system. In order to carefully account for statistical 
weights, a Boltzmann weight equal to w is considered only if the bend happens: 

1. At the mth level and at least one polymer bond at the level is inside the square. 

2. At the (m + l)th or (m + 2)th level, and both polymer bonds at the level are inside 
the polymer. 

A weight w p = w a is considered for any configuration in which two bonds are parallel 
to each other within the same square. We can, furthermore, distinguish configurations in 
which two bonds are parallel to each other from configurations where three consecutive bonds 
form a hairpin configuration. As explained above, whenever this configuration is present, 
an additional weight w\ Y = w b is introduced. 

It is not important to know along which of the two lower bonds in the (m + l)th or 
(m + 2)th square does the polymer chain enter into the mth square. In fact, even if the 
polymer undergoes a bend while moving from the higher level square to the mth level square 
the corresponding weight is already taken into account into the partial partition function of 
the higher level site. 

A weight w c = w c is considered for any configuration in which a filler particle and a 
monomer belonging to a polymer chain are nearest neighbors. In order to consider two 
particles in contact we need them to be nearest neighbors on the lattice. In other words, 
it must be possible to reach one particle from the other by moving one single step on the 
lattice. For example, the base and the peak sites inside one plaquette are not nearest 
neighbors because there is no single lattice bond connecting them. The same is true for the 
two middle sites of the square. The base site or the peak site is a nearest neighbor of the 
two middle sites. 
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C. Recursion Relations 

Following the usual procedure explained in Ref . [| , the next step consists in constructing 
;he recursion relations (RRs) for the partial partition functions. Contrary to Refs. Q and 
3], we start right away with constructing the RRs for a particular 2-cycle solution, suitable 
for describing the ground state (at absolute zero), in which the base site of the plaquette 
under investigation has label m, but we make a distinction between the middle sites and 
the peak site so that the middle sites have label m + 1 and the top site has the label m + 2. 
There is no harm in doing this as the 1-cycle solution is a special case of this scheme as 
noted in the previous paper of the series, and as we will further clarify below. In other 
words, the scheme will produce a 1-cycle solution that will describe the disordered phase at 
high temperatures. Therefore, it is a more general solution that a 1-cycle solution. 

The recursion relation for Z m (l), the partial partition function of the mth branch of 
the Husimi lattice given that the mth level site is in the I state, is written as a function 
of the Z m+ i(a) and Z m+2 (/3), the partial partition functions describing the contributions 
from the upper branches and summing up all the contributions coming from all the possible 
configurations that are compatible with the state at the mth level site. 

Z m (l) = w 2 w h w p Z m+1 (0)[Z m+1 (L)Z m+2 (R) + Z m+1 (R)Z m+2 (L)} 

+ wZ m+1 (R)Z m+1 (L)(Z m+2 {Y) + wlZ m+2 (S)). (3) 

Considering the case in which the mth level site is in the O state, the partial partition 
function Z m (0) for the O state can be written as: 

Z m (0) = Z 2 m+l {\)Z rn+2 {\) + 2w 2 c Z m+1 (l)Z m+2 (l)Z m+1 (S) + w 2 c Z 2 m+1 (l)Z m+2 (S) 

+ 2w 2 Z m+1 (l)Z m+2 (S)Z m+1 (S) + wlZ 2 n+1 {$)Z m+2 {\) + w 2 Z^ +1 (S)Z m+2 (S) 
+ [Z m+1 (T) + w 2 Z m+1 (S)] [Z m+1 (L)Z m+2 (R) + Z m+ i(R)Z m+2 (L)] 
+ wZ m+1 (R)Z m+1 (L)Z m+2 (0). (4) 
When the mth level site is in the L state the partial partition function can be written as: 
Z m (L) = [Z m+1 (R) + wZ m+1 (L)}{Z m+1 (T)Z m+2 (T) + w 2 Z m+1 (S)Z m+2 (I) 

+ w 2 Z m+1 (l)Z m+2 (S) + w 2 Z m+1 (S)Z m+2 (S) + w 2 w p w h Z m+1 (0)Z m+2 (0)} 
+ w p [wZl l+1 (L)Z m+2 (R) + Z 2 m+l {R)Z m+2 {Vj\ 

+ [Z m+2 (R) + wZ m+2 (L)] wZ m+1 {Y)Z m+1 {0). (5) 
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The relation for the R state is easily obtained from Z m (L) by the interchange L < — > R: 

Z m (R) = [Z m+1 (L) + wZ m+1 (R)}{Z m+1 (l)Z m+2 (T) + w 2 Z m+1 (S)Z m+2 (T) 

+ w 2 Z m+ i(S) ^ m _|_2(S) + w 2 WpWi 1 Z m+ i( K 0) Z m+ 2(0)} 
+ w v [wZ^ +1 {R)Z m+2 (L) + Z 2 m+l (L)Z m+2 {R)] 

+ [Z m+2 (L) + wZ m+2 (R)] wZ m+1 (l)Z m+1 (0). (6) 

Finally, the recursion relation for Z m (S), the partition function of the mth branch of the 
Husimi lattice given that the mth level site is in the S state, is given by: 

Z m (S) = r][w^Zl l+1 (T)Z m+2 (T) + 2w 2 Z m+1 (T)Z m+2 (l)Z m+1 (S) + w^Z^ +1 (l)Z m+2 (S) 
+ 2w; Z m +l(l)Z m + 2 (S)Z rn+ i(S) + W c Z m+1 (S)Z m+2 (I) + (S)Z m+2 (S) 
+ \w 2 c Z m+1 (\) + w 2 Z m+1 (S)] [Z m+1 (L)Z m+2 (R) + Z m+ i(R)Z m+2 (Vj] 
+ ww 2 c Z m+1 (R)Z m+1 (L)Z m+2 (0)]. (7) 

It is possible to write analogous relations for Z m+ \{a) by properly substituting m — > m+ 1, 
m + 1 — >• m + 2 and m + 2 — > m + 3. 

We introduce B m = (Z m (L) + Z m (R)) and the following ratios between partial partition 
functions at the mth level site of the lattice: 

x m {a) = Z m {a) I B m , a = I,0,L, and S 

x m (R) = 1 - x m (L). (8) 

We divide the partial partition functions by B m since we are interested in the limit in which 
at low temperature the entire lattice is covered by the polymer. So, B m is never zero, and 
the ratios remain finite at every temperature. 

As we climb down any of the branches of the Husimi lattice, a site can be classified as 
a simultaneous base site (of the upper square) and peak site (of the lower square). These 
sites differ in their levels (that is, the lattice generation they belong to) by two. Thus, it 
should come as no surprise that as one moves from a level that is infinitely far away from 
the origin towards the origin itself, the recursion relations can approach 2-cycle fix-point 
(FP) solutions, x 2m (a) — > ao, x 2m+ i(a) — > a±, etc., where a = I, O, L,R or S, regardless of 
the initial values of the ratios at infinity. The indexes and 1 refer to even and odd levels, 
respectively. At the fixed point, the values of a and ai become independent of m. Thus, 
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these fix-point solutions of the recursion relations describe the behavior in the interior of 
the Husimi tree. We express the fix point in this case as follows, where m is taken to be an 
even integer including 0: 

x rn (T) = x m+2 (l) = z , x m (0) = x m+2 (0) = oq, x m (L) = x 2m+2 (L) = Z , 
x m (R) = x m+2 (R) = 1 - x m (L) = 1 - Z , x m (S) = x m+2 (S) = s 
x m+ i(I) = x m+3 (I) = n, x m+1 (0) = x m+3 (0) = oi, x m+ i(L) = x m+3 (L) = ii, 
x m +i(R) = x m+3 (R) = 1 - x m+1 (L) = 1 - k, x m+1 {S) = x m+3 (S) = si. (9) 

The system of equations for the ratios can be obtained from the system of equations 
for the partial partition functions and it is easy to show that they can be written in the 
following form: 

i Q = w(l- k) hi + w 2 w p w h o 1 [k(l - l ) + Z (l - h)] , (10) 
°oQ = + 2w 2 Sii 1 i + w 2 s il + 2w 2 SiiiS + w^sfio + w 2 sls 

+ h [h(l - Z ) + / (1 - h)} + w(l - h)ho , (11) 
l Q = (1 - h + wli)[iii + w c 2 si« + w^iiSo + w^sis^ + w 2 w p w k OiO Q ] 

+ woiii{\ — lo + wl ) + w p (l - k) 2 l + ww v l\(l - Z ), 
s oQ = vWl^o + 2wlsiiiiQ + w c s il + 2w 2 s 1 i 1 s + w 2 s 2 i + s 2 s 

+ K 2 zi + w 2 Sl ] [h(l - Iq) + Z (l - Zi)] + un/7 c 2 (l - Zi)Z l0o , (12) 
iiQ' = w(l- Z ) Z z'i + w 2 w p u7 h oo [i a (l - *i) + Zi(l - l Q )] , (13) 
OiQ' = i&i + 2w 2 s i ii + w 2 c s h i 2 Q + 2w 2 c s i Si + w^o^i + w c s l s i 

+ 2 [/o(l - h) +h(l- k)] + w(l - Z )Z 0Ol , (14) 
ZiQ' = (1 - Z + wl )[i ii + w 2 s ii + Wc^oSi + w 2 s s h + w 2 w p w h o oi] 

+ w;o io(l - h + wli) + w p (l - / ) 2 /i + ww p ll(l - h), (15) 
SiQ' = vWl^lh + 2w 2 s io^i + w t s i^l + 2w 2 s ^oSi + w c 2 s^i + SqSi 

+ K 2 i + w 2 s } [Z (l - Zx) + Zx(l - Z )] + w c 2 (l - Zo)Zoo x . (16) 

where Q' is obtained from Q by exchanging and 1 subscripts and Q is given by: 

Q — (1 + w)[Mo + ^c-Mo + W^1 S + ^c s l s + W 2 W p W h OiO ) 

+ ujoiii(l + w) + io p [(l - /i) 2 Z + w^/ 2 (l - fo)]- (17) 
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The polynomials Q and Q' relate B m with B m+ i for even and odd levels as follows: 

B m = Bf n+1 Q, (m even), B m = B 3 m+l Q, (m odd). (18) 

At low temperatures, the only solutions that we have been able to find is the 2- 
cycle solution in which the ratios at even levels are different from those at odd lev- 
els. At high temperatures, this difference between even and odd levels vanishes so that 
%m( a ) — x m+ i(a) = x m+ 2(a), a = I, O, L, R or S, and the 2-cycle solution reduces to a 
1-cycle solution, as noted above. 



D. Free Energy 

In order to determine which phase is the stable one at any temperature, we have to 
find the free energy corresponding to the grand canonical partition function |2 of all the 
possible phases of the system as a function of w. The free energy per site, which is related 
to the osmotic pressure, can be easily calculated from the expressions for the total partition 
function Z at the (m = 0)th, (m = l)th and (m = 2)th levels Q]. 

The total partition function of the system Zq can be written by considering the two 
(m = 0)th branches Cq meeting at the origin, the (m = 0)th level site. For this, we need to 
consider all the possible configurations in the two branches. This is done by considering all 
the configurations that the system can assume in the two squares that meet at the origin of 
the lattice. There are two configurations which corresponds to states in which a bond gets 
to the origin but does not cross it, each of them contributing 

Zo(I)Z (0) (19) 

to the total partition function. This term reflects the fact that such configurations look like 
on I state as looked from one side of the origin and as an O state as looked at from the 
opposite side. Two other configurations correspond to cases in which a bond passes through 
the origin and makes a bend at the origin, each contributing 

(l/™)Z , g (L)Z , g (R), (20) 

where the factor (l/u>) is needed in order not to take into account the Boltzmann weight 
for the polymer bend at the origin twice and the subscript "g" refers to the gauche part of 
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the partition function for L and R states. In fact, it is possible to separate a "gauche" and 
a "trans" contribution to the partial partition functions for the R and L states at any level. 
The "gauche" portion is the one corresponding to configurations such that there is a bending 
at the mth level site, while the "trans" portion is the one corresponding to configurations in 
which the two bonds coming out of the mth level site that we are considering are straight. 
It is easily seen that: 

(L) — Z m+ i(R){Z m+ i(l) Z m+2 (1) + w c Z m+ i (S)Z m+2 (I) 

+ w 2 Z m+1 (T)Z m+2 (S) + w 2 Z m+ i(S)Z m+2 {S) + w 2 w p w h Z m+1 (0)Z m+2 (0)} 

+ w p Z 2 m+1 (R)Z m+2 (L) + Z m+2 (R)wZ m+1 (l)Z m+1 (0). (21) 

+ w 2 Z m+1 (l)Z m+2 (S) + w 2 Z m+1 (S)Z m+2 (S) + w 2 w p w h Z m+1 (0)Z m+2 (0)} 

+ w p wZ 2 n+1 (L)Z m+2 (R) + wZ m+2 (L)wZ m+1 (l)Z m+1 (0). (22) 

Two configurations correspond to cases in which a bond passes through the origin but 
does not make any bend at the origin, one contributing 

^t(L) (23) 

and the other one contributing 

Z 2 o, t (R) (24) 

depending on the position of the bonds at the origin. Finally, one configuration corresponds 
to the presence of a solvent molecule at the origin, with a contribution: 

(1/V)Z (S) 2 . (25) 

Also in this case, the factor (l/rj) is necessary in order not to not to take into account the 
activity for the filler particle at the origin twice. It is then possible to write: 

Z = 2Zo(I)Z (0) + (2/™)Z , g (L)Z , g (R) + Z 2 t (L) + Z 2 t (R) + Z 2 (S)/ V 

= B 2 Q 2 , (26) 

where B = Z (L) + Z (R), and where we have introduced the polynomial 

2 

Q 2 = 2i o + -lo, s r , g + /Jo + r t,o + (27) 
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where r ,t, ?"o,g^o,t^o,g are the trans and gauche portions of r and I, respectively. 

The total partition function can be used to obtain the thermodynamics of the system. It 
is clear that Z is the total partition function of the system obtained by joining two branches 
Cq together at the origin. Now, let us imagine taking away from the lattice the two squares 
that meet at the origin. This leaves behind six different branches C\. We connect these 
branches to form three smaller but identical lattices, the partition function of each of which 
is denoted by Z\\ the latter can be written in a form that is identical to that of equation for 
Zq, except that the index of each partial partition function is replaced by 1. 

The difference between the free energy of the entire lattice and that of the three reduced 
lattices is just the free energy corresponding to a pair of squares. This pair of squares 
contains one complete site, that is the center of the lattice, and six additional "half" sites 
where each of the remaining C\ branches has a base. So, these two squares correspond to 
a total of four lattice sites. Following Gujrati P, |2| we can write the adimensional free 
energy per site without the conventional minus sign as 

* = i>n(§). (28, 

It is possible to write 

Zi = BfQ 2l (29) 

similar to the way we had expressed Z = B^Q 2 earlier, and where B m and Q 2 are defined 
above. Since B = BfQ, see (jTSjl . the free energy per site can be written as 

u = ±h l (Q-Y (30) 



2 \Q2, 

In order to determine which phase is the stable one at some temperature, we must find 
the free energy of all the possible phases of the system as a function of w. The usual 
thermodynamic potential per site, which is the osmotic pressure, can be obtained from uj 
through: 

F = -Too. (31) 



E. Phase Diagram 



The typical phase diagram that is obtained is shown in Figure 01 The values of the 
parameters of the problem for this particular case are shown in the figure caption. 

1(3 



There are three phases that are captured by this calculation, a crystalline phase and two 
liquid like phases. The crystalline phase is captured by the 2-cycle recursion fix-point while 
the two liquid phases are obtained using the 1-cycle recursion fix-point. The crystalline phase 
(CR) is stable at low temperature. If the system is in the CR phase and it is heated up it 
melts into the liquid phase at the melting temperature Tm- The liquid phase that is stable 
above such transition is called the equilibrium liquid (EL) in order to distinguish it from 
a second liquid phase that is always metastable (and that we label as ML, the metastable 
liquid). 

If we start from high temperature in the EL phase and crystallization is avoided at the 
melting temperature, the system finds itself in the supercooled liquid phase. This phase 
merges continuously into the ML phase at a temperature Tmc that has many characteristics 
of the critical temperature predicted by the mode-coupling theory 22|. Using the model of 
polymer that we have introduced above, it turns out that that the density of filler particles 
is zero in both the CR and ML phases. These two phases have an ordered structure such 
that the filler is pushed out of the system and all the sites of the lattice are occupied by 
a monomer belonging to the polymer. The filler density is non zero in the EL phase. The 
filler density as a function of the temperature for the system whose free energy is shown in 
Figure El is shown in Figure 0] It is very easy to observe in the figure how the filler density 
goes to zero as the temperature approaches T M c- 



F. Percolation 



It is then interesting to study the percolation process in the EL phase keeping in mind 
that this phase is the stable one above Tm and it is metastable below Tm and above Tmc- 
The effect of the change in the chemical potential of the filler particles on the phase diagram 
is minimal. The values of Tm and Tmc are changing very weakly with the value of the 
chemical potential so that we focus on one particular representative set of values for a, b 
and c. 

In order to calculate the filler percolation we proceed as follows (see Ref. [lj). Following 
the Gujrati approach js^ ]. we introduce the probability R m < 1 that a site occupied at 
the mth generation is connected to a finite cluster of occupied sites at higher generations. 
Then, Z m R m denotes the contribution to the partial partition function C m due to all those 
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configurations in which the site at the mth generation is connected to a finite cluster in C m . 

If we divide Z m R m by Z m , we obtain a recursion relation for R m . At the fixed point the 

value of R m does not depend on m anymore, and we denote it by R. For this particular 

problem, the equation for R is: 

w 2 i 3 + 2Rw 2 i 2 s + w^i 2 s + 3R 2 w 2 is 2 + R 3 s 3 + 2lr(w 2 i + Rw 2 s) + ww 2 olr 
w 2 i 3 + 2w 2 i 2 s + wf:i 2 s + 3w 2 is 2 + s 3 + 2lr(w 2 i + w 2 s) + ww 2 olr 

We can then obtain the percolation probability from the value of R by means of 

p = l-R 2 (33) 

where we consider percolation to have occurred only if the central site of the lattice is 
connected to an infinite cluster on either or both sides of the lattice, as explained previously 

DD 

We have studied the dependence of the strength of the percolation process on three 
variables: the chemical potential of the filler particle, the degree of semiflexibility of the 
polymer and the strength of the interaction of filler and polymer. The results are presented 
in Figures El to |H1 In all the figures, the filled dots correspond to temperature intervals where 
the disordered phase in which the fillers are embedded is thermally stable while the empty 
symbols correspond to intervals in which it is metastable, that is between Tm and Tmo 
Figure El shows the behavior of the percolation probability of the filler particles as a function 
of temperature for different values of the chemical potential of the filler. An increase in the 
magnitude of the chemical potential makes the percolation less likely to occur. 

Figure El shows that the temperature range at which the percolation of the filler particles 
occurs becomes narrower and narrower as /x becomes more and more negative. In particular, 
for values of fi larger in magnitude than 0.5, no percolation occurs in the system at any 
temperature. Above some critical value of the chemical potential, the filler particles do not 
percolate anymore. From an energetic point of view, if the chemical potential of the filler 
particles becomes too large and negative it becomes highly unfavorable for the system to 
have filler particles in the equilibrium configuration and the density of the filler drops to 
values that are too low for the percolation to occur. The tailoring of the chemical nature 
of the filler, as a function of the polymer matrix to be used, seems to be of the uttermost 
importance. 

Figure El demonstrates that the strength of the percolation process can still be very strong 
in the metastable phase. In the temperature intervals corresponding to the empty symbols 
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in the figures, the system is not percolating in the ordered stable phase while it is percolating 
in the disordered metastable phase. We have observed something very similar in a parallel 
study [2I when discussing the percolation in the case of a system made of particles of different 
sizes (but with no polymer matrix present in the system). In that case also, the metastable 
phase presents a stronger percolation than the ordered phase. These observations suggest 
how it might be useful to prepare the system in a metastable phase in order to obtain a 
stronger percolation. From an experimental point of view, a metastable phase can be very 
similar to a stable one in many aspects. As long as the lifetime of the metastable phase 
is very long with respect to the observation time or the expected lifetime of an application 
produced using such a metastable phase, it might be useful to have a system in a metastable 
phase in order to increase the strength of a percolation process. We have also studied the 
percolation process as a function of the semiflexibility of the polymer chain. As we have 



argued elsewhere 



24j and above, the parameter a = e p gives a measure of the stiffness of 



the polymer chains. The value a = corresponds to the least flexible polymer chain, and 
the original Flory's model. As the value of a increases, the polymer chain becomes more 
and more flexible. This is easy to understand if we remember that increasing a is equivalent 
to decreasing e while keeping e p constant. 

Figure |H1 shows the value of the percolation probability as a function of temperature 
for different values of a. An increase in the stiffness of the polymer matrix makes the 
percolation of the filler particles more and more unfavorable. In particular, while for a < 0.1 
the system is not percolating anymore at the lowest temperature at which the metastable 
phase exists, for larger values of this parameter the system is still percolating when it gets to 
Tmc- Even larger values of a would just move the percolation threshold to higher and higher 
temperatures. The stiffness of the polymer matrix apparently opposes the percolation of the 
filler particles. While there is no percolation in the ordered crystalline phase, the increased 
stiffness makes the percolation more and more difficult to occur in the disorder phase as well, 
in both the stable and metastable regions. We are not aware of any experimental systematic 
study of this kind of dependence: one might argue that a decreased stiffness could help the 
percolation process by creating some preferential channels for the filler particles to percolate. 
However, this would be misleading because the dependence on (repulsive) interaction e c must 
also be taken into consideration in addition to the chemical potential. The chosen values of 
the two conspire in a such a way to the above result. 
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Finally, we have studied the dependence of the percolation process on the strength of 
the interaction between polymer matrix and fillers. In the model, the parameter c = e c /e 
is associated with the excess interaction between the filler particles and the monomers be- 
longing to the polymer chain. An increase in c corresponds to an increase in the repulsion 
between the filler particles and the polymer matrix. The results summarized in figures [7| and 
IHJshow that the percolation process becomes more favorable as the repulsion between filler 
and matrix increases. This is particularly clear in Figure |H1 where it is possible to observe 
how a change of c from 0.04 to 0.08 induces a change in the percolation threshold from 0.25 
to 0.245. When the repulsion between matrix and filler is larger, it takes a lower density 
of filler in the system to obtain a percolating cluster. This should not come as a surprise 
as increasing the repulsion will eventually bring about a phase separation between fillers 
and the polymer matrix, in which the the solvent phase will be percolating. This kind of 
phase separation is not desirable for many application. The above-mentioned increase in 
percolation tendency is a precursor of this eventual phase separation in the model. 

As the repulsion becomes weaker, the filler particles are better distributed in the matrix 
and they tend to form smaller and smaller clusters since they have almost no preference 
between sitting next to another filler particle or a polymer unit. Thus, from an experimental 
point of view, it seems necessary to introduce filler particles with quite a strong repulsion 
towards the matrix in order to increase the probability of having percolation in the system 
if a stronger tendency to percolation is desired. 

Thus, it is clear once again the tailoring of the filler to the matrix appears to be funda- 
mental to obtain optimal properties. 

III. PERCOLATION OF DIFFERENT-SIZED FILLER PARTICLES IN THE 
PRESENCE OF A POLYMER MATRIX 

In order to study the effect of size and shape disparity on percolation of fillers in a 
polymer matrix, we have considered the addition of two different particles to the single-site 
fillers whose percolation has been described in the previous section. We have studied the 
percolation process in the following systems: 

• single-site particles and square particles in the presence of a polymer matrix 
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• single-site particles and star particles (occupying 5 lattice sites) in the presence of a 
polymer matrix 

This choice allows us to draw some conclusion about the effects of the size and shape 
of the particles on the percolation properties. In the following, we will always refer to the 
single-site species as the "solvent" particles even when they are not the majority component. 

A. States of a site 

In order to solve the problem recursively, we have to define the different states a site can 
be in. In the case of square and single-site particles, every site of the lattice can be in one of 
seven states. If we look at an mth level site, this site can be occupied by either a monomer 
that belongs to a polymer chain and that can be in one of the states I, O, L or R, introduced 
in the previous section, or by a solvent molecule, in the S state, or by the edge of a square, 
either in the A or in the B state, according to the position of the square with respect to the 
site The site under investigation is in an A state if a vertex of the square is at the site 
and the square lies above the mth level. The site is in the B state if the vertex of the square 
is at the mth level site but the square lies below the site. 

In the case of star and single-site particles, every site of the lattice can be in one of eight 
states. If we look at an mth level site, this site can be occupied by either a monomer that 
belongs to a polymer chain and that can be in one of the states I, O, L or R, introduced 
above, or by a solvent molecule, in the S state, or by a star, either in the A or in the B or 
in the C state Q]. As said above, we consider stars that occupy five sites: a central one, 
named as the "core" and four endpoints that are all nearest neighbors of the core but (on 
the square and Husimi lattices) are not nearest neighbors of each other. The A, B and C, 
states are defined as follows |2|: the site under investigation is in an A state if one end point 
of the star is at the site but the core of the star is above the mth level. The site is in the 
B state if an end point of the star is at the mth level site but its core is below the site. 
Finally, the site is in the C state if the core of the star is at the mth level site that we are 
considering. 

We now consider the two cases separately. 
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B. Squares and solvent molecules 



The total partition function of the system made of solvent molecules and square particles 
on the Husimi lattice is 



The number of solvent molecules N$ is controlled by the solvent activity 77 = where 
fj, is the chemical potential of the solvent particle. The number of square molecules N sq is 
controlled by the activity 1% = u> _/ib , where /Xb is the chemical potential associated with 
one side of the square particle. There are three different interactions between nearest neigh- 
bor contact pairs N c , N e between solvent particles and monomers, square particles and 
monomers and solvent particles and square particles, respectively. The corresponding Boltz- 
mann weights are w c = exp (— /3e c ) = w c , = exp (—/fed) = w d and w e = exp (— /fe c ) = w e , 
respectively, where the e's are the corresponding excess energies. The parameters that de- 
scribe the nature of the polymer are the same ones introduced in the previous section. Since 
we are interested in the case in which the system is a pure polymer phase in the ground 
state at T = 0, we always consider negative values of /1 and /ib- 

The procedure to obtain the free energy and the percolation probability is the same as 
explained in Ref. 3]. The complete recursion relations are very cumbersome but straight 
forward to derive by following the procedure shown explicitly in the first two papers 
What one needs to do is to consider all possibilities for a plaquette to determine their con- 
tribution to any partial partition function. As the problems becomes complex, the number 
of possibilities grow rapidly and patience is required to ensure that all terms are accounted 
for; several tens of them are present in some of the equations. We give the final recursion 
relations in the appendix for the benefit of a reader. 

In all the cases that we have studied, both in the presence of square and star particles, 
the qualitative features of the phase diagram have not changed. In all of the results, we 
always have three polymer phases: a crystalline one stable at low temperature, an equi- 
librium liquid phase stable at high temperature and a second metastable liquid phase that 
is never the stable one but that is still very important since the metastable extension of 
the equilibrium liquid below the melting temperature of the system, what we call the su- 
percooled liquid phase, merges into the metastable phase at T MC . The free energy versus 




(34) 
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temperature diagram for all these systems looks qualitatively the same as that shown in 
Figure H3 with minimal quantitative corrections to the values of Tuc and Tu so that, once 
again, we just focus on one representative set of parameters. The density of both solvent 
and square particles in both the crystalline and the metastable phase is identically zero and 
the percolation can only take place in the liquid phase, both in its equilibrium branch and 
in the supercooled extension below the melting temperature. 

When we are looking at the interactions between the different particles and their perco- 
lation properties, there are two possible approaches: 

1. We can consider the solvent particles and the square particles as different species with 
a nonzero excess energy and consider the separate percolation of the two kinds of 
particles in the presence of the other ones. 

2. We can consider the solvent and square particles to be of the same chemical nature 
and consider the total percolation probability when particles of different sizes and 
same chemical nature are present in the system. Here the cluster is formed by both 
particles. One can also study the percolation of each kind of particles separately. 

The first question that can be asked is the following: does a system containing just square 
particles exhibit any percolation? In order to study this we have set the activity for the 
solvent molecules equal to zero and studied the system of square particles embedded in a 
polymer matrix. 

The results of this calculation are the following: regardless of the value of their chemical 
potential, the square particles alone never percolate at any temperature even if their chemical 
potential becomes extremely small in magnitude. We think this might be related to the 
nature of the polymer we have considered here and of course is also related to our choice 
to always consider negative values of \i and fit, in order to have a pure polymer phase in 
the ground state at T = 0. Having an infinite molecular weight polymer chain makes it 
almost impossible for these bulky particles to form an infinitely large percolating network. 
By considering a system of polymer chains with a finite degree of polymerization, we might 
have percolation of these larger particles as well. One might also think about considering 



a less restrictive condition of connectivity 



calculation carried on during this research 



or percolation. Here, as well as in any other 



ELQ 



we have considered two particles to be part 
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of a cluster only if they are nearest neighbors of each other on the lattice. It is possible to 
define the condition for percolation in less strict terms: for example it is possible to consider 
percolation of second nearest neighbors or even of particles farther apart with respect to 
each other on the lattice. However, we will not do that here. 

The first problem to study is the problem of the percolation of the complete system in 
which we consider the solvent particles and the square particles to be of the same chemical 
nature and just different in size. To some extent, we can think about this system as a system 
of particles with a bimodal distribution of sizes. In order to consider this problem, we set the 
excess energy for the solvent-square interaction equal to zero and we consider, out of all the 
possible configurations of the system at the mth level, those configurations corresponding 
to a percolating cluster of solvent and square particles. 

Figure 01 shows the percolation probability for the complete system as a function of the 
temperature. The values of the parameters of the system are shown in the figure caption. 
One immediately notices how the temperature interval in which the system is percolating 
is much wider than in the case of the percolation of the single-site filler particles studied in 
the previous paper If we consider the dependence of the percolation probability on the 
total density of filler, the corresponding results are shown in Figure ITU1 

The value of the percolation threshold is much lower than in the case of the single-site 
fillers alone, for the same value of the chemical potential (per site) for the particles. Thus, 
the results show that size disparity lowers the percolation threshold and makes percolation 
more likely to occur, even in the presence of a polymer matrix. 

It is also possible to study the effect of the presence of the larger particles on the per- 
colation of the smaller ones as compared to the case in which only the smaller particles are 
present in the system. 

Figures E] and ^] show the percolation probability for the smaller filler particles as a 
function of temperature and filler density, respectively, for three different cases: in one case, 
only the smaller filler particles are present; in the second case both particles represent but 
they do not have any interaction; and in the third case both the particles are present and 
they have a repulsive interaction between themselves. 

It is clear from the results how the presence of the larger particles makes it easier for 
the smaller particles to percolate, lowering their percolation threshold. The effect of the 
interactions between the different particles seems weaker since we do not see a significant 
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difference between the curves corresponding to the cases in which the two kinds of particles 
interact and do not interact between themselves. Even if the effect is small, though, it is 
clear how the presence of repulsive interaction between the particles slightly enhances the 
percolation process, as seen in Figure IT21 where the percolation threshold is slightly lower in 
the case in which repulsive interactions are present between the particles. Figure IT31 shows 
how much the percolation probability is increased by the presence of the larger particles. 
The filler density is, in both cases, the total density associated with all the filler particles 
present in the system. 

C. Stars and solvent molecules 

For a system made of solvent particles, star particles and polymer, the total partition 
function is identical to that for squares except that N sq is replaced by the number of star 
molecules N st This number is controlled by the activity ??b = u>~ Mb , where /ib is the chemical 
potential associated with an arm, which also happens to be a bond, of the star particle. All 
other terms in have the same meaning as above, except that stars replace squares. Since 
we are interested, as before, in the case in which the system is a pure polymer in the ground 
state at T = 0, we always consider negative values of /i and /ib- 

We give the final recursion relations in the appendix for the benefit of the reader. 

As in the case of the problem of square particles, the first question that we have asked is 
the following: does a system made of just star particles in the presence of a polymer matrix 
exhibit any percolation? In order to study this problem, once again we have set the activity 
for the solvent molecules equal to zero and studied the system of star particles embedded in 
a polymer matrix. 

Exactly like in the case of square particles, we have obtained that regardless of the value 
of their chemical potential, the star particles alone never percolate at any temperature. As 
explained in the previous subsection, we think this might be related to the nature of the 
polymer we have considered here. 

As in the previous case, we have studied the problem of the percolation of the complete 
system in which we consider the solvent particles and the star particles to be of the same 
chemical nature and just different in size. In order to study this problem, we set the excess 
energy for the solvent-star interaction equal to zero and we consider, out of all the possible 
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configurations of the system at the mth level, those configurations corresponding to a per- 
colating cluster of solvent and star particles. Figure ITU shows the percolation probability for 
the complete system of solvent, stars, and polymer, as a function of the temperature. The 
values of the parameters for this system are shown in the figure caption. We immediately 
notice how the temperature interval in which the system is percolating is much wider than 
in the case of the percolation of the single-site filler particles but at the same time much 
narrower than that of the system of squares and solvent molecules studied in the previous 
section and shown in Figure |H1 The results for the dependence of the percolation probability 
on the total density of filler are shown in Figure The value of the percolation thresh- 
old is much lower than in the case of the single-site fillers alone, for the same value of the 
chemical potential (per site) for the particles. Thus, the results show how the presence of 
size disparity, in this system as in the system containing the square particles, lowers the 
percolation threshold and makes percolation more likely to occur. 

It is also possible to study the effect of the presence of the larger particles on the per- 
colation of the smaller ones as compared to the case in which only the smaller particles are 
present in the system, as done before for the case of the square particles. Figures EI] and El 
show the percolation probability for the smaller filler particles as a function of temperature 
and filler density, respectively, for three different cases: in one case only the smaller filler 
particles are present, in the second case both particles are present but they do not have any 
interaction, and in the third case both the particles are present and they have a repulsive 
interaction between themselves. 

A careful observation of Figure El and its comparison with Figure E3 shows a very 
interesting result: while the addition of the square particles has a minimal effect on the 
percolation properties of the smaller particles and induces a slight decrease in the percolation 
threshold for the smaller particles, the addition of star particles increases the percolation 
threshold for the smaller particles making percolation more unlikely to occur. This difference 
may be related to the fact that these two particle shapes that have been used, squares and 
stars, have a different size but the same number of nearest neighbors interaction that they 
can create with other species, eight in both cases on a Husimi lattice. It is clear form the 
results how the presence of the larger particles makes it easier for the smaller particles to 
percolate, lowering their percolation threshold. As we have done before, we have to point out 
that at the present stage it is more proper to state that this might not be the only important 
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factor affecting the percolation threshold and that a host of properties might be involved. 
The effect of the interactions between the different particles seems not very important since 
we do not see a significant difference between the curves corresponding to the cases in which 
the two sets of particles interact and do not interact between themselves. Even if the effect 
is small, though, it is clear how the presence of repulsive interactions between the particles 
slightly enhances the percolation process, as seen clearly in Figure IT7I where the percolation 
threshold is slightly lower in the case in which repulsive interactions are present between the 
particles. Figure ITH1 shows how much the percolation probability is increased by the presence 
of the larger particles. The filler density in both cases is the total density associated with 
both the solvent particles and the larger ones (either squares or stars). 

It is very useful to compare the percolation probability for the two composite systems that 
we have described in these two last sections. Figure shows the percolation probability 
as a function of the total filler density for the two systems that we have investigated. It is 
clear how the percolation threshold is lower in the case of the system containing the square 
particles as compared to the case of the system containing the star particles. It is possible to 
argue that the star particles are larger than the square particles since the former occupy five 
sites of the lattice (the core and the four end-points) while the square particles occupy only 
four sites. As observed in experiments, an increase in polydispersity and in the size disparity 
of the filler particles is usually associated with a decrease in the percolation threshold and 
our model calculation is able to capture this feature. 



IV. CONCLUSIONS 

We have developed a lattice approach that allows us to study the effects of size and 
shape of particles and their energetics on the percolation properties of a composite system. 

The approaches that have been used so far in the literature have a wide range of lim- 
itations that range from the use of random mixing approximation to the need to have 
incompressible systems to the requirement for the monomers and voids to have the same 
size. The most important limitation of all existing theories is represented by the necessity to 
go through some averaging process on scales that are much larger than the size of nanoscopic 
fillers found in many of the most advanced composites that are currently under investigation. 
Molecular dynamics and Monte-Carlo simulations can be used in order to obtain a better 
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understanding of the thermodynamics of the problem; however, both these techniques are 
extremely time consuming and require finite size scaling that is not obvious in the case of 
complex systems. Also, simulations usually require an unreasonably large amount of free 
volume and cannot be carried out in the nearly incompressible limit. 

Our approach allows us to describe how percolation is affected by the presence of size 
and shape disparity between the species that are present in a system. This approach is very 
powerful since the calculations are straightforward and exact. The solution that we obtain 
represents an approximate theory on the original lattice. This kind of calculation represents 
a new method to deal with nanoscales where all the interactions can be taken into account 
and where the nanoscopic details of composite systems can be easily taken into account. 

Our results are able to reproduce most of the experimental findings that have been 
observed in the literature. In particular: 

• the percolation threshold decreases as the size disparity between different particles 
present in the system increases. We have proved this by studying systems made of 
single-site, square shaped or star shaped filler particles; 

• the percolation threshold decreases with the aspect ratio of the particles present in 
the system. We have obtained this results through the comparison of the behavior of 
systems containing either the star or the square particles; 

• the stiffness of the polymer matrix as well as the nature of the interactions between 
fillers and matrix affect the percolation properties of the system; 

• in many different systems it is possible to observe the presence of a metastable phase 
that is percolating very strongly while the stable phase has a very weak percolating 
cluster (a much lower value of p) or even no percolating cluster at all. This observation 
suggests that it might be useful to prepare a system in a metastable state in some cases 
in which a percolating network of particles is needed in order to enhance one or more 
physical properties of the composite system. 

There are still numerous aspects of this problem that are worth investigating. The effects 
produced by branching and finite loops have not been investigated so far to the best of our 
knowledge. One aspect that will be very interesting to study is the effect of the finite chain 
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length of the polymer matrix on the percolation of the filler particles. When the polymer 
chains are short, either linear or branched, they can give rise to gelation, the formation of an 
infinite network of polymer chains. The gelation of the polymer chains and the percolation 
of the filler particles are two independent processes but it is easy to understand how these 
two processes might strongly affect one another. Our polymer model has been recently 
implemented in order to describe the presence of finite polymer chains 2^. The current 
model describes a polydisperse polymer system with an equilibrium degree of polymerization 
but the model can be easily extended in order to be able to describe a monodisperse polymer 
system. In the case of the infinitely long polymers studied here, the larger particles never 
percolate in the absence of the small ones when they are inserted in the polymer matrix. 
We feel that this could be related to the nature of the polymer chains used and that the 
inclusion of finite length polymers might in principle produce new and exciting results. 
Another natural extension of this model will involve the study of the percolation process in 
the presence of even larger particles even though it has to be kept in mind that an increase 
in the size of the particles will lead to an increase of the complications to be faced in the 
description of the system. 
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FIGURE CAPTIONS 
FIG. 1: Portion of an infinite lattice known as Husimi tree. 

FIG. 2: Possible configurations of the site at the mth level of the tree for a system containing 
polymer chains and monomeric solvent molecules: I state, O state, L state, R state and S state. 

FIG. 3: Free energy as a function of temperature for the polymer system with monomeric filler 
particles. The curves for the crystalline phase, the equilibrium liquid and the metastable liquid 
phases are shown for a = 0.20, fi = —0.50, 6 = 0.07. 

FIG. 4: Solvent density as a function of temperature for the polymer system with monomeric 
filler particles. The curves for the crystalline phase and the equilibrium liquid are shown for 
a = 0.20, n = -0.50,6 = 0.07. 

FIG. 5: Dependence of the percolation probability for monomeric filler particles in the presence of 
a polymer matrix on the temperature of the system for different values of the chemical potential 
of the filler particles. For all the curves: a = 0.1, 6 = 0, c = 0.07. 

FIG. 6: Dependence of the percolation probability for monomeric filler particles in the presence 
of a polymer matrix on the temperature of the system for different values of the flexibility of the 
polymer chains. For all the curves: /i = —0.2, b = 0, c = 0.07. 

FIG. 7: Dependence of the percolation probability for monomeric filler particles in the presence of 
a polymer matrix on the temperature of the system for different values of the interaction energy 
between filler particles and polymer matrix. For all the curves: [i = —0.2, a = 0.2, 6 = 0. 

FIG. 8: Dependence of the percolation probability for monomeric filler particles in the presence 
of a polymer matrix on the volume fraction of filler for different values of the interaction energy 
between filler particles and polymer matrix. For all the curves: /i = —0.2, a = 0.2, 6 = 0. 

FIG. 9: Dependence of the percolation probability on the temperature for the system containing 
monomeric filler particles along with square particles. Here fx = —0.2, fj,^ = —0.2, a = 0.2, 6 = 
0,c = (i = 0.07,e = 0. 
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FIG. 10: Dependence of the percolation probability on the total density of filler for the system 
containing monomeric filler particles along with square particles. Here \x = —0.2, \x\> = —0.2, a = 
0.2, = 0, c = (i = 0.07, e = 0. 

FIG. 11: Effect of the addition of a second filler species and of the interactions between fillers 
on the percolation probability of a composite system as a function of temperature: squares and 
solvent. 

FIG. 12: Effect of the addition of a second filler species and of the interactions between fillers on 
the percolation probability of a composite system as a function of the filler density: squares and 
solvent. 

FIG. 13: Percolation probability as a function of the filler density for the system containing 
monomeric fillers and square particles compared with the percolation probability of the system 
containing only the smaller particles. 

FIG. 14: Dependence of the percolation probability on the temperature for the system containing 
monomeric filler particles along with star particles. Here /i = —0.2, ^ = —0.2, a = 0.2, b = 0, c = 
d = 0.07, e = 0. 

FIG. 15: Dependence of the percolation probability on the total density of filler for the system 
containing monomeric filler particles along with star particles. Here [i = —0.2, fi^ = —0.2, a = 
0.2, b = 0, c = d = 0.07, e = 0. 

FIG. 16: Effect of the addition of a second filler species and of the interactions between fillers on 
the percolation probability of a composite system as a function of temperature: stars and solvent. 

FIG. 17: Effect of the addition of a second filler species and of the interactions between fillers 
on the percolation probability of a composite system as a function of the filler density: stars and 
solvent. 

FIG. 18: Percolation probability as a function of the filler density for the system containing 
monomeric fillers and star particles compared with the percolation probability of the system con- 
taining only the smaller particles. 

FIG. 19: Percolation in systems made of partigles of different shape: effect of the difference in 
particles sizes and shapes on the percolation probability of the system as a function of the filler 
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